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Abstract. Cmpact general formulae for the energy, momentum and angular 
momentum radiated by confined systems of charges and currents are presented in 
terms of their multipole Cartesian tensors. Besides the usual electric and magnetic 
multipoles, a family of toroid (anapole) electromagnetic tensors, as well as some mean 
squared radii contributions are standed out through a method given in previous works 
of the authors. 



1. Introduction 

The purpose of the present paper is to obtain general formulae for the energy, momentum 
and angular momentum loss in the case of a bounded charge and current distribution, 
in the form of multipolar expansions, using only fully symmetric and trace-free (STF) 
electric , magnetic and anapole (toroidal) multipole tensors, together with some mean 
squared radii. The gauge invariant technique of reducing the usual Cartesian electric 
and magnetic multipole moments to STF ones , presented in [2] , is used for this end. 
In the section 2 the multipolar expansions of the scalar and vector potentials, $ and 
A, are presented in terms of electric and magnetic Cartesian multipole moments. In 
section 3 the coresponding expansions of the electric and magnetic fields, E and B, are 
given. The relevant terms for the purposes of this paper are only those that fall off as 
fast as 1/r and 1/r 2 as r — > oo given in equations (fTTj) and (fl"2"|) . In section 4 the total 
radiation intensity (energy loss ) is calculated and given by equation (flUj) in the form of 
an expansion in terms of STF multipole moments . The same calculation is performed 
in section 5 and 6 respectively for the total momentum loss and angular momentum loss. 
In section 7, the approximation's criteria are discussed and, finally, one presents explicit 
expressions for the radiated energy, momentum and angular momentum including all 
the electric and magnetic multipoles up to order three (octopoles). These results are 
compared with some incomplete expressions existing in literature. 

In appendix A one can find the definition formulae for the reduced electric and 
magnetic multipole moments . 
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2. General Formulae 

We shall accept in the following a compromise regarding the unit systems used for writing 
Maxwell's equations and the different definitions and equations for electromagnetic 
theory. Let the Maxwell's equations be written as 

a a at a at 

VB = 0, V-E = —p (f) 

where the positive parameters Eq, Hq, a are constants depending on the system of 
units. Considering the consequences of Maxwell equations, from physical reasons, only 
two independent constants will be needed because of the constraint 

— = c 2 (2) 

c standing for the speed of light in vacuum. 

If the Maxwell equations are written in the form above, than the fields E, B and 
the potentials <£>, A are related by the equations 

1 dA 

B = V x A, E = - V$ — (3) 

a dt w 

while the retarded potentials are given by 

V V 

One can write the equations and formulae above in the International unit system 
by taking a = 1 and the IS values for e , fiQ. The same relations can be written in 
Heaviside-Lorentz unit system by taking e = 1, /z = 1, ex = c in cgs units while in 
Gauss unit system e = 1/47T, p,o = 47r, a = c . 

The goal of this compromise is to compare easily the results from the present paper 
with similar results in different issues written either in IS or in Gauss or Heaviside- 
Lorentz units. 

Let us consider a charge p(r,t) and a current j(r,t) distributions having supports 
included in a finite domain T>. Choosing the origin O of the Cartesian axes in V, the 
retarded vector and scalar potentials at a point outside V, r = x^ei^ will be given by 
the multipolar expansions 

-or Ft V (~ 1 ) n ~ 1 A a Mii-in-i.kjto) -, 
^-i^ijk^j , u %\ ■ ■ ■ u % n -\\. J 

n>l "" ' 

+ -eiE^^^---^A Pil "^ liit0 \ t = t- r -, (4) 
a n\ r c 
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teo4(r , t) = E t2L v» . « = E *. • • • *. [5^1 • 

The electric and magnetic moments are defined as 

P<»)(t) = / £>(^)d 3 £ : Pix..*. = / & • ••C,/^./!<l :i C (5) 

M w (f) = 7 / r x j(U)d 3 £ : M,,., n = / & . . x jkd 3 £ 

m + l)a j n+1 qj 

v v 

It was shown in [2] that we can perform such transformations of the multipole tensors 
pW—.pW M W^M W , (6) 

(where P , M are STF tensors) so that, if A and $ are obtained from equations 
(jlj after the substitutions (jBj), the correspondence (A, $) — > (A, $) is a gauge 
t r ansf or mat ion . 

In Appendix A are presented the principal results of the reducing procedure. In the 
following we will suppose that all the multipole moments are represented by reduced 
STF tensors. 

3. Expansions of the fields E and B 

Let us write the multipole expansion of B = V x A. Use of equation (jlj) and formulae 
V x (V x a) = V ■ (V • a) - A a and 

[A - ± * ]«^& = 0, for r + 0, 
cr at 1 r 

leads to the expansion 

B(r,t) = ^t]p v . [V n^ ito)]] 

<±/l n > 1 ft. f 

_ E (zlTl v»- . [ifi w W ] + / £ -v x £ tip- v- . [ip'-Vo)], 

4ttc 2 ^ n! L r v ;J 47ra n! L r v ;J ' 

to = t- T -. (7) 
For i2 we get 

In order to obtain the behavior of the fields E and B at large distances, we need the 
formula for the partial derivatives of a retarded arbitrary solution of the wave equation. 
The following formula can be easily obtained: 

dh--- di n [-f{to)\ = -——^i • • • v, 



, 1 (-1)" 



r v /J r c n 1 " dt" 

jn-l x/^ 

0j 



fan-l tn in-2-in-linJ 



+ 0(4) (9) 
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where 

D n = £> n _i + n, D = 0, v = -. (10) 

r 

In equation Q and in the following we understand by Au x ...j n \ the sum over all the 
permutations of the symbols i ? that give distinct terms. The use of the formula Q in 
equations (J7J) and (jSJ), and also of the well known properties of STF tensors, allows the 
ordering of the contributions in the expansions above as a 1/r power series. 

+ £i E i K«K • - A,!-"- 1 • - ^ x K- . I*']} + • • .(ID 

with the notation / >n = d n /(to)/dt n and 

Let us write £ = E^ + E^ +0 (1 / Y 3 ) where -E^ 1 -* and -E*- 2 -* are respectively, proportional 
with 1/r and 1/r 2 and, similarly, B = B (1) + B (2) + 0(l/r 3 ). The parts E {1) and B (1) 
of the field are purely transverse fields, satisfying the properties 

v ■ = 0, u ■ B« = 0; EM = -B« x , £ |£ (1) )| 2 = -|B«| 2 . (13) 

a /io 

4. Radiation Intensity 

The Poynting vector will be written in terms of E^ and fiW by considering the 
equations (|T3j): 

S = - (E x B) = e \E^\ 2 cv + 0(h = -\B^\ 2 cu + O(l) 
/i H /i H 

The total radiated power may be written as the limit of the integral on the sphere of 

radius r as r — > oo: 

2 

J= hm J r 2 v SdO(i/) = | rV-SdUH = ^~ / l^ (1) | 2 d^(^) 

= — <r 2 |B( 1 )| 2 > (14) 
where 

We obtain from equation (|TT|) that 

r 2 |B«| 2 = (^) 2 £ -^(-^(^-MtU^.M^) 

v 4vr y n ^ 1 n!m! I c n+m+2V ' n+1/v ' m+i; 
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■ «!2») • • m!:' +1 ) + -^-^ . mI:' +1 ) . \v x (»— . pi:' +I )i 

The following formula is used to perform the average operations, PQ 

f 0, n = 2k+l, } . s 

< y i • • • ^ >= i a. ., n - 2 fc jfe-0 12 (16) 

I (2fc + l)!! 1*1*2 • • • "Jn-l^nj! "' — Z,ft J ft — U, J., Z,, . . . I 

If A*™), B*" 1 ) are STF tensors, one can show that 

< (^.AW).(/.B< m ') >= T^Yjn ( A(n) - B(m) ) 

For the calculation of the averaged square of the vector .B^ one can apply the following 
formulae obtained from equation (|17jl : 

< K.M (n) )(^.M M ) >= -^%-(M (n) .M W ), 
v A ; (2n + l)!! v ; ' 

V -M " x (^.P W )]) = ((^^M^W^.P^), 
= < . . . i n - 1 v i v h . . . Uj m _ x > M il ... in _ li P il „. im _ lfc = 0. (18) 

By introducing the expansion (fTTj) of -B in equation (jl4j) . and using the relations 
(fTHJ) . one obtains 



a 2 n + 1 



'M (n) M (n) ) + ^(p {n) p (n) ) 



(19) 



47T£ c 3 ^ nn!(2n + l)!!c 2n 
5. Recoil Force 

By considering the momentum current density tensor as 

%i = l(e Q E 2 + —B 2 )5 l3 - (soEiEj + —B^j), (20) 
I Ho Ho 

then the recoil force is given by 

F R = - lim * (u.T)dS = - lim r 2 / ^,dQ(i/) e,- 

= - lim r 2 / -( £o £ 2 + — B 2 )v&n(v) = - — ( r 2 \B^\ 2 v). (21) 
r ^°° •/ 2 // Ho 

Let us introduce the expansion (f!3j) of l-B^I 2 in equation (j2~Tj) and consider the relations 

n]S m,n-l fijin) ^j(n-l)^ (w + 1) !<? m , w +l j^^K 

(2n + l)!! 1 n+1 n j+ (2n + 3)!! 1 n+1 ' '"+ 2 j ' 
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(n - l)!5 mn _i ~(n) rj(n-l) N , ra!5 mjn+ i ~(n) r^(n+l)v 

= (2n-l)H ( M .M-i- M ,n )+(^TI)!! (M ^- M ^ ) ' 

(K^.m^V^x (,™- 1 .p (m V> = K-vw-O 

~ ~ (ra-l)! ~ 

x e i £ gj'fcMil...i n _iqPjX---jm-lfc = — ^77, -\- \\\\ *1— *n-19^*l— im-lk^nm 

(n — 1)! /r^W ~(n). „ , \ 

e^ fc (M V ; .n-l.P l J )J m , (22) 



~~ (2n+l)!! 3 V Jjk"nm, 
with the obvious notation, 

(M l \n-l.P { ) jk = Mi X) ...i n -ijPii...i n -ik- 

After observing that the terms that give nonzero contributions are only those containing 
a 6 factor, one gets 
F = Afo 1 
* 47r^c 2 "(2n + l)!! 

f™ + l/rs( n ) n + 2 ,~(n) ri(n+i). 

X {^ (M " + " M " ) + C 3(2, i + 3)H(, i + l)! (M ^- M ^ ) 

2 -e« ,.(M W „-lP W ] e(" + l) ( pW p(»-D, 



acnin J crn! 



a 2 c(2n + 3)!!(n + 1)! 1 ' n+1 ' n+2 j J 1 dj 

Since M^ -* = 0, and the total electric charge in T> is constant (dP^ -* /dt = 0), then the 
following formula can be obtained after performing some changes of summation indices: 

p __V0 V 1 f n + 2 ~(n) ~(n+l) C 2 ~(n) -(„+!) 

27rc»4i C *» \(n + l)!(2n + 3)!! [ n+1 n+2 + a 2 1 ' n+1 " ' n+2 J 

4 n!n(2nVl)!! e ^ ( ^-" " (24) 
6. Angular Momentum Loss 

The total angular momentum lost per unit time by a radiating system is given by the 
flux of the radiated electromagnetic angular momentum through the spherical surface of 
radius r —>■ oo. The 2nd-order antisymmetric tensor associated with the electromagnetic 
angular momentum density is defined by its Cartesian components m^- = XiQj — Xjgi 
where the electromagnetic momentum density vector is g = (eo/oe)(E x B). The 
angular momentum density current is associated with the 3rd-order tensor jJkij — 
XiTjk — Xj%k with Tjk given by equation (|2U|). By introducing the angular momentum 
density pseudovector L with the components Li = (l/2)eyjj,mjfe, and the 2nd-order 
angular momentum current tensor /i^ = (l/2)e K y//jS, we obtain that 



(2) 



e (u.E)(u x E)i + —iv.BMv x B) 
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r— >oo 



and further that 
V dt /rad 

= 47rr 3 ^(i/.B (2) )(^ x +e (f -£ {2) )(f x 

By using equations (fTU|l - (fT2Jl . one obtains 



r 



^.B( 2 ) = /rac/i 47r^ n + 



and 



i/ n .M 



n! C" 



E^ 



1 



r „ x B W = J^L 

4:ira c' l n\ 



a 



-v x (v 



n-l 



M 



(n) , 
■ n+lJ 



K.py + > + (*/ 



n-l 



j(n) v 



(25) 

(26) 
(27) 

28) 



ri/ x E W = 



Aira c n n 



E x (^.p^i) + k-m;: +1 ). - (^.m;: +1 [(29) 



(n) 



We than have the following substitution rules: 



M 



in) 



;(«) 



4p w 

After using equations (J22|) and 

-i)> 



((i/\M (n) )(i/ m-1 .P (m) NN - ra!5m -™ +1 



(2n + l)!! v ' n ' n+2 ' 



(30) 



(31) 



and the substitution rules (jSOJ)? we obtain 
47rr 3 ((^.B( 2 ))(^xB«)> = ^£ 

^ n>l 



n + 1 



;(n) 



n!(2n + l)!!c 2n + 



T e^ jfc (M in |n- 1|M 



n + 1 



(n- l)!(2n + l)!!^ 2 ™" 1 



n + 2 



n!(2n + 3)!!ac 2n+1 



v^,n • P,n+2 , 



(32) 



and 



4vrr 3 ((i/. J B( 2 ))^x = ^£ 

^ n>l 



e i e ijk{r^ n \n - 1|H n+ iJ jA . 



n!(2n + l)!!a 4 c : 



2n-3 



n + 1 



:p ( ?.M ( r 1K 



n + 2 



(n-l)!(2n + l)!!a 3 c 2 "- 3V ' n ' n y n!(2n + S)!!^ 2 ™" 1 
By introducing equations (jHZj) and (jHHjl in equation (j2HJ), we obtain] 



(33) 



'd£r\ _ //o 

, dt /rad~ 47r nSl 
n + 1 

~ ac 2n-i( n - i)!(2n + 1)!! 



n+1 



c 2n+1 n!(2n+ 1)!! 



— ( n ) 

ej£ iifc (M n .n - l.M 



n + 2 



ac 2n+1 n\(2n + 3)\\ 



(n) . 

n+l)jk 



p(n+l), 



a 2 c 2n-l n! ( 2n+ 1)1! 

n + 1 

ac 2n-l( n _ l)!(2n+ 1)!! 



,pM jyj(n-l), 

v. , n ,n j 



v,,l r ,n - IVI ,n+2 ) 



ac 2n+1 n!(2n + 3)!! 



(34) 
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ie algebraic operations, we obtain 

J_ / n+1 ^ 
- C ^\ n!(2n + l)!! ei%fcX 

.n-l.P ( ;M +(mJ; ) .ti-1.M^ 1 ) + (35) 

/ jk \ / j'fe 

+ ,r +2 „, [p ( :i 1} • m ( ;i - m ( £? ■ p ( ;ii +^: ] ■ - p ( n n) ■ m^I } . 

an\(2n + 3)!! L i n ,n+2 j | 



dL 




dt 


47TC 


X 


cr V ' 



7. Discussion and conclusions 

The procedure, given in the present paper, for calculating the energy, momentum and 
angular momentum radiated by arbitrary sources is characterized by the mathematical 
simplicity. In our calculations only simple algebraic manipulation is necessary. The 
results expressed in different coordinate systems, for example spherical, may be obtained 
by adequate transformations. Complete and correct multipole analysis is given in 
some fundamental works as (HI El • I n these works, the multipole expansion formulae 
are introduced using the canonical basis of the solutions of the scalar wave equation 
of Helmholtz. Although the results obtained by this formalism represent an exact 
description of the multipolar expansion applied to the electromagnetic field of arbitrary 
sources, the formulae are rather cumbersome. Due to the complexity of these formulae, 
the possibility to obtain erroneous results exists even in the case one calculates the 
contributions of the low order multipoles to the radiated field. This is the reason why 
we believe it is necessary to give corrections for the radiation formulae obtained in [Bl Ej 
in addition to the corrections given in these papers to the results from [H1E1- 

By applying the formulae for the radiated quantities, some care is necessary for 
obtaining a correct grouping of the different multipolar terms in a given approximation. 
For understanding this problem it is sufficient to refer to the harmonic time variation 
of the charges and currents. Let us consider 

p(r,t) = p (r)e-^, j(r,t) = joMe"^. 

By considering the magnetic multipole tensor M(t), the n-order time derivative may be 
written as 

d^ M H = (-iw) w n r ( )d 3 ^ ^ (36) 

dt n (n + l)aJ y J v ; 

where d is the linear dimension of the domain T>. By introducing the wave length of the 
emited radiation, A = 2ttc/uj, we may write 

I£MM~(!V. (37) 

The parameters u and d are related, as order of magnitude, to the charges velocities, 
ujd ~ v, such that ud/c ~ v/c. In what concerns the electric multipole's contributions, 
we notice that the contribution of the time derivative of the n-order electric multipole 



Multipolar radiation and the gauge invariant reduction of multipole tensors 



9 



is of the same order of magnitude as the contribution of the (n — l)th order magnetic 
multipole because 

d P h ... in (t) = J j- V(& • • • &„)d 3 £ = J • • • C„. ,/;. ;<l :i C 



dt 



v v 



Consequently, as is pointed out also in |5], in the case of the radiation intensity, in 
equations (fT9|) . (}24|) and (}35j) . a consistent expansion should take into account, alongside 
with the magnetic multipoles up to a given order /i, the electric multipoles up to the 
order e = 

For a comparison with some results from the literature, for example [HJ U\, we 
give below the results of the expansions of the total radiated power, recoil force and 
angular momentum for (//, e) = (4, 5) pointing out the criteria for grouping the different 
multipolar terms for a given approximation. 

In order to understand this problem, we give for the case of the radiated power a 
detailed description of the procedure used. We denote by X^ ,£ ) the radiated power 
obtained by the recursive procedure begining from the order /i for the magnetic 
multipoles and from the order e for the electric ones. In the case of the harmonic 
time variation, we stand out only the parameters d, c and u. In the case (/x, e) = (4, 5), 
considering the c dependence of the coefficients given by equations (|A.20|) . and using 
equations (|A.22|) - (|A.26j) . we can write 

~W U 2 d 2 UJ A d A ~( 2 ) , U 2 d 3 

P ~ a 10 + a n —^- + Q i2-^ + • • • > p ~ a 20 d + a 21 —^- + . . . , 

~(3) uj 2 d 4 ~( 4 ) . ~(5) 

P ~ a 30 d 2 + —g- + . . . , P ~ a 40 d 3 + . . . , P ~ a 50 d 4 + ...; 
M ~a 1Q d + a n — — + ..., M ~ a' 2Q d z + a 21 — — + ..., 

~(3) ~(4) 

M ~a' 30 (i 3 + ..., M ~aV 4 + ... (38) 

where and a'^ are functions of time proportional with sin tut or cos cut, and the 
magnitude order of a common factor is not considered. We point out that in equations 
(JHBj) a contribution to the absolute order of magnitude of the final results gives the 
current density j ~ v ~ ud. This factor present in all the expressions from equations 
(JHEJ), associated with 1/c factors from equations (|T9|) . (j2lj) and (}35|) . gives a factor 
(d/X) 2 ~ {v/c) 2 in the expressions of the corresponding expansions. This common factor 
is irrelevant four our considerations. If we consider equations ()38|) in the multipolar 
expansion (|19|) . we can write in the case (/j, e) = (4,5) the following result for the 
relative orders of magnitude of the different terms 

1^ ~A Q + A 2 ( 2 + A 4 C 4 + A 6 ( G + As( 8 , C = d/X. (39) 

From equations (|A.15|) and (|A.16|) . one can see that X (5 ' 6) can be written in the form 
of a similar expansion in which A , A 2 , and A 4 are unchanged, and only Aq, A 8 ,...will 
be different. Consequently, for a consistent approximation of I using the multipole 
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expansions up to (fi,e) = (4,5) one should keep only the terms of the order (d/X) n for 
n < 4. 

Similar considerations may be applied for the evaluation of the recoil force and 
radiated angular momentum. From the results given in jH] and (7j, it seems that the 
authors claim to give the first 3 terms of the d/X expansion of the radiated power, 
the first two terms of the expansion of the recoil force and again the first 3 terms of 
the expansion of the angular momentum loss. Our evaluations give for these cases, 
very different results. It seems that, in the above mentioned papers, the approximation 
criteria are not applied consistently. The comprising of the factor 1/c in the definitions 
of different parameters as, for example, for the toroid dipole, seems to be one of the 
sources of errors. We point out the benefit of using a system-free expression of Maxwell's 
equations, as in the present paper. So, the factors c included in different definitions in 
Gauss system of units is represented here by a which is no counting to the evaluation 
of approximation order. This fact was observed also in [3] in the case of the energy loss. 

We give bellow our results for the above mentioned approximations. 

In the case of the radiated power, it is necessary to calculate Z^ 4,5 ) and retain only 
the powers of ( up to four. One obtains 



1 



4tte c 3 



2.. 2 
-P 



3c2 P ' * 4 20c 2 



10c 



1 ... --(2) 

V . T (1) + 



2 --(3) 

V 



945c 4 



1 -.(2) ...(2) 2a 2 .. 9 2 -2 4 . ^ 



3c 2 



2 ,2 

3c 1 



■■■■(3) 4a 2 .... a 2 

V H rm. u H t 

3c 4 ^ 20c 4 



... (2) 

M • M 



(2) 



+ 0(C 6 



(40) 



The recoil force is given by 



11 

AtteqF = — - -p.V 

C 
1 -(2) ... 



■ ■■(2) 2a 



-m x p 



2a 
3c 2 



H xp — rax t 



Of .. - (2) 



(2) 



4 -(2) 



5c 2 yt " ' (1) ' " ' ' ' ' 315c 2 
and the radiated angular momentum 



V 



•(3) 



a 



30c 2 



('-,:, jk MqjVqk -0(C 



AlTEr 



dL 

~dt 



! 2 ... 2 2a 2 a 
-p x p H -px t H -t X p -rh x rh H -rh 

V 3c 2 ^ 3c 2 y 3c 2 5c 2 



,p(2) 



a . -(2) 



^e i£i:jk P qj V qk ~p x T (2) - — tx t 



2a 2 ... 2 - .. •>•>-' 



3c 4 



3c 4 
a 

5c 4 
1 



2a' . ■■■■ a ... 
— -mx. 11 H j 11 .V [2) j 

3c 4 ^ 5c 4 ^ 5c 4 



a .. -(2) 
m. T (1) 



5c^ 



5c 4 J 



a . :r( 2 ) a ... ~(2) 



5c 4 



(41) 



.V.(2) 



a •■ ... (2) 

2a ...(3) ... (2) 
V . M 



^i^-ijkPqj qk j^g^4 ^i^'ijk ~^qjl~^qk 



a- 



10c 



-eiEijkAiqj Aiqk 



+ 



105c 4 
2 

315c 4 



^i^ijk *P qpj'P 'qpk 



2a ...(3) 
105?^ 

+ o(c 6 ). 



...(2) 2a —(3) .. , , 



2 « •- (3) ~f2l 



105c 4 



(42) 
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The following notations was introduced in the last three equations: 

p, m, /x, t for 7>« T$, and T (2) = ej^, (43) 

For comparison, we give below the results from [HI EI with our notation but written in 
Gauss unit system. So, the result from is | 

2 ... 1 - ,2 2 . .. |2 1 . ... 1 «: ... 1 «: 

^dt = - - 2 t I + — |m| + — ( -- T (1)ij )( Vi, -~ 2 T m ) 

2 

+ 945c7 ( MijkMijk + MijkMijk )■ (44) 

In this expression are present terms as AdijkAiijk and ~T of order six in £ and 
are missed terms of the order four as VijkVijk for example. This expression may not to 
be considered as an expansion up to the order six in £ beacause in this case many terms 
must be added. 

The result for X from [7] is given by 

2 ., 2 ., 1 2 -2 4 

1 1 , , 

VijVij +7^ MijMij ■ (45) 



20c 7 3 J 20c 7 

—2 

In equation ()45jh the terms i and m. /x represent corrections of order four but many 
terms of the same order are missing compared with equation (|40|). Similar conclusions 
may be formulated regarding the results for recoil force and angular momentum. 



Appendix A. Gauge invariant reduction of multipole Cartesian tensors 

The transformations © involve a sequence of operations meant to obtain the symmetric 
and traceless part of the tensors implied. Let an n-rank tensor of magnetic type 
i.e. symmetric in the first n — 1 indices and verifying the property \-i 1 ...i k „ 1 ji k+1 ...i n -ij — 
0, k = 1 . . . n — 1. The the symmetric part of this tensor is given by 

1 1 n ~ * 

L(sym)ii...i n = "~[Lii...t n '-i rt i 2 ...ii + • • • + Lji...j n — \-i\...i n J]] ^ixinl^h—in-i q \}-~ ^] 

71 U A=l 

where J\f-( x > is the component with the index i\ suppressed. 

The operator Af defines a correspondence between lJ n ) and a tensor of rank (n — 1): 

l>> — Af[l»] . [AT[L^]] . . = M,..^, = £ i „_ lps U 1 ..,„_ 2ps (A.l) 

L Jli...I„_i 

which is fully symmetric in the first n — 2 indices and the contractions of the last index 
with all the previous indices give null results. So, the tensor A/"[L^] is of the type 
ft/|(n-i)_ Particularly, 

-l) k n 



2 lM(n)] = / x jd% 

(n + 1 a J 



(n + l)a 

I For comparison with our formulae one takes ^tteq — y 1, a 



V 



Multipolar radiation and the gauge invariant reduction of multipole tensors 12 

A ^2fc+l[ M (n)] = (-1) " J ^2*^n-2fc-l x ^ x ^) d 3^ = 0,1,2. . . (A.2) 



V 



Let a fully symmetric tensor and the detracer operator X introduced in 0j. This 
operator acts on a totally symmetric tensor so that T'fS^- 1 ] is a fully symmetric 
and traceless tensor of rank n. The detracer theorem states that [3] § 

rrrsWll -T (-lH2n-l-2m)!! (n:m) 

where [n/2] denotes the integer part of [n/2], Afa...^} is the sum over all permutations 
of the symbols z 9 which give distinct terms, and S^ m "^| in denotes the components of 
the (n — 2m)th-order tensor obtained from S^ n ^ by contracting m pairs of symbols i. 
It is useful to introduce here another operator A by the relationship 

T[SW]1 = V, n -% i2 [A[SW]] (A.4) 

where A[S^] defines a fully symmetric tensor of rank n — 2. From this definition 
together with the theorem (jA.3J) . we obtain 

A- ■ rsWl-'y (-l) m ^-l-2(m + l)]!! (n:m+1) 
^H-t n -2 P J- (m + l)(2n — 1)!! 2 ' ' ' l2m - ll2mi 2 m +i-i n -2}- 

In the following, for simplifying the notation, all arguments of the operator A should 
be considered as a symmetrized tensor i.e. A[T (n ' 1 ] = A[T^J for any tensor T^. The 
same applies to the operator T: T[T( n )] = T[T^]. 

The following four transformation properties of the multipole tensors and potentials 
may be used for establishing the results from 0. 

I. Let the transformation of the nth-order magnetic tensor: 

M<"> - Mg : M [L)il .. An = M h ... in - ijXi^L-i ff [L (B) (*o)]- (A.6) 

71 A=l 

Let us substitute in the expansion of the potential A the tensor M^ n ^ by M/?% obtaining 



(L) 

"0 (-I)™" ' 



1 (L)J 4tt n\n 3 1 

- (^ 1 M,.,„_ li + . . . 8 jin _ 1 N il ... in _ 2i ))[\S n \t Q )) 



^ n— 1 

- E ^iife^i A fc 9 A/'/ 1 A) in _ 1(? [L (n) (i(:o)] 



A=l 



A + ej c^c^ . . . <9j n _ 1 [(5ii 1 A / i 2 ...j n _ 1 j + . . . o~a n _ 1 N'i 1 ...i n _2j / 



A + Vtf(r,t) + ^ ( - 1) " ( "" 1) V-' 
47T n'crn 



(f )] 

r 



(A.7) 



Here the equation [A — (l/c 2 )9 2 /9t 2 ][/(t )/^] = 0, r 7^ is considered. The function 
\I/ is a solution of the homogeneous wave equation for r / and the corresponding 

§ In this equation, the definition of the symmetric and traceless part of the tensor S^ 71 * differs from 
that used in 0] by a factor l/(2n — 1)!!. 
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expression is irrelevant. Let the transformation 

p(n-l) ^ p'(n-l) = p(n-l) + 0l ( n ) < /^[|_W] j fll ( n ) = _. °LV—^1. (A.8) 



c 2 n 2 

Introducing the transformed potentials produced by the substitution p( n_1 ) — > p'( n - 1 ) ; 
we obtain 

A[mW->mW p(n-i)_>p'(n-i)] = A + V*, $[pc«-i)-*p'(»-i)] = $ - — . 

So, the transformation (|A.6|) produces changes in the potentials which, up to a gauge 
transformation, are compensated by the transformation (|A.8|) . 
II. Let the transformation of the nth-order electric tensor: 

n-l 



(A.9) 



A=l 



We obtain 
A[p(")-»pW] = A + 



Mo {-l) n -\n - 1) 



47T 



V x V 



n! n 



m-2 



: A/"[L (n) ] 



$[p(»)-»pW] = $. 



The change of the vector potential is compensated by the transformation 



M (n_1) — > M (n - 1} + a 2 (n)AT[L (n) ], a 2 l 



77 



77 — 1 



— a x (n). 



(A.10) 



(A.11) 



III. Let the transformation of the magnetic vector of rank n: 
M (n) _ M W . M(s)ii ^_. n = M n .., n - £ ^A^SWft,)] 

where is a fully symmetric tensor. The alteration of the vectorial potential is 
eliminated by the transformation 

77-2 



M(™- 2 ) — ► M (n - 2) + b(n) A[S (n) ], 6(77) 

IV. The transformation 
p(») — > p. 



2c 2 n 



(S) 



P(S , )«i...j„ — Pij...in E ^1112^3 -in [** ] 

Z>(<) 

produces the following changes of the potentials: 



(A.12) 
(A.13) 



A[pw^pW] = A 



/'o 



47r n! 



0(0 



7 a«o (-i) n ~v-2)(n-i) vn - 3 

47r 2n\c 2 



- A [S (n) l 

r 



with as satisfying the homogeneous wave equation and 



^ r , m ^ Mo (— 1)" 1 n(n — 9 

1 (S)J 4tt 2c%! 

Let the transformation 

p("-2) > p»(«-2) = p(n-2) + 6 ( n ) A [S(t Q )] 



-A[S (n) ] 

r 



(A. 14) 

with 6(n) given by equation (|A.12|) . The effect of the transformation (jA.14|) on the 
potential A is the compensation of the extra-gauge term . So A[p< n )^p&!,p(' l_2 )-»-p''( n_a )] = 
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A + V\P' but it is easy to see that the modification of the scalar potential $ 
produced by the transformation (jA.14|) together with the modification due to the 
transformation ()A.13j) give $[pM^p(™)[s( n '],p( 71 - 2 )^p"("- 2 >] = §-d^' /dt the total effect of 
the transformations (jA.13|) and (jA.14|) being a gauge transformation of the potentials. 

4. Let the gauge invariant process of reducing the multipole tensors begin from the 
order n = e in the case of the electric tensors and from n = fi for the magnetic ones, 
and go downward up to lowest value, n — 1. Generally, we may choose e > \i as will be 
seen in the following. In are given the formulae for the results P and M : 



P 



I [(e-»»)/2] V2k 

?><»>+ r E 4 n) ^ k [p {n+2k) ] 



[(M-n-l)/2] ,2k+l k , 

+ E ^rrE^A^-^lM^ 1 ^] , (A.15) 
fc=o ar z=o 1 

, . f[(n-n)/2] ,2k k ) 

M (n) = M {n) + r\ E Taf E tf« 5 A'.A/™ -2 ' [M (n+2fc >] I (A.16) 

[ fc=i i=0 J 

where 

4 n) = n^+2o, 



B kf = HKn + 2q) J[ ai (n + 1 + 2k - 2h) JJ a 2 (n + 2A; - 2s), (A.17) 

o=l /i=0 s=0 



?= 

i fc-i-i fc-j-i 

CjM = [] 6( n + 2g) JJ ai (n + 2k - 2h) J[ a 2 (n - 1 + 2k - 2s) (A.18) 

q=l h=0 s=0 

and 

, . a n — 1 . . 1 n — 1 c 2 . N , , . n — 2 

Oi(n) = — - — a 2 {n) = — = - a^n), b(n) = — — . (A.19) 

c z n z a n z a z 2c z n 

odi 



c ib a n a, zc ib 

By introducing equations (jA.19|) in equations ()A.17|) and (jA.18|) . we obtain 

. iwi 1 n 



H ~ 2 k c 2k n + 2k' 

|(») _ 



-l) k ~ l+l u n{n + 2l)\ 



M 2 l c 2k + z (n + 2k + \){n + 2k + 1)!' 

2<c 2fc (n + 2fc)(n + 2fc)!' 1 ' ; 

7>( n ) and M {n) stand for the "static" expressions of the reduced multipole tensors: 

V (n \t) = T[PW] = i~\T7 / p(r,t)r 2n+1 V n -d 3 

(2n — 1)!! J r 



* W <*> = T|M( " ,, = a(n + 1 K2n- 1)!! g 1 / r2 " + '^'-' f ) * 
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In formulae above one should consider that rifc=z Fk = Hi L < /.In the case e = /i + 1, 
we may write 



,(«) 



l(e-n)/2] 



k j2fc-l 



•l) fe d 



-T 



(n) 



k=l 



T (n) 



i)V fc r 



c 2fc (J^2fe-1 (fc)' 
fc-1 



(A.21) 



4 n) A fc [P (n+2fc} ] + £ B^AW 21 " 21 ' 1 ^ 2 ^} 

1=0 

where besides the usual electric and magnetic multipole moments, a third multipole 
family, the toroid moments and, generally, mean-square radii of various orders, are 
involved [UJ [7j . 

For applications, we we give the results of the gauge invariant reduction of the 
electric and magnetic multipole tensors, begining the reduction procedure from the rank 
e = 5 for the electric tensors, and from // = 4 for the magnetic ones. From equation 
(IA.21|) we have for the electric multipole tensors, 



(1) i t (D 1 ? CD 

' c 2 '(1) + c 4 '(2), 



[T 



(1)1 

(i)k 



\J {1) ] 



T {AWA[P (3) ] + 2^[M< 2 >]} = ^ / [(* ■ i)6 - 2^]d 3 e, 



,( 2 ) 



7?(2) 



1 -j-(2) 
72 ' (1)' 



[T 



(2)1 
(1) life 

1 

T2 



c 2 T ^4 2) A[P 



,( 4 ) n 



,(3) 



P (3) _ Ljf] 



[T 



(3)i 
(1)1 ikl 

1 

60 



c 2 ri4 3) A[p (5) ] + 4 3 ^[M (4) ] 



ikl 



V 



,(4) 



p ^ _ p(4) p( 5 ) = £>(5) 

The magnetic multipole tensors are given by 



M 



(i) 

i 

(2) 
ik 

(3) 



M - A^ (1) + T |Ci; W 2 [M w ] + C#A[M 
MS' = + T feX[M (4) ] + C gA[M (4) 



,( 3 )! 



,(1) 



(3), 



1 

mi + -^fii 



ik 



(A.22) 



(A.23) 



(A.24) 
(A.25) 



(A.26) 
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where 



m = 




Ua J 



42a 



1 



(A.27) 



v 



V 
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